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ON THE SUBRIEMANNIAN GEOMETRY OF CONTACT ANOSOV FLOWS 


SLOBODAN N. SIMIC 


Abstract. We investigate certain natural connections between subriemannian geometry and hy¬ 
perbolic dynamical systems. In particular, we study dynamically defined horizontal distributions 
which split into two integrable ones and ask: how is the energy of a subriemannian geodesic shared 
between its projections onto the integrable summands? We show that if the horizontal distribution 
is the sum of the strong stable and strong unstable distributions of a special type of a contact 
Anosov flow in three dimensions, then for any short enough subriemannian geodesic connecting 
points on the same orbit of the Anosov flow, the energy of the geodesic is shared equally between 
its projections onto the stable and unstable bundles. The proof relies on a connection between 
the geodesic equations and the harmonic oscillator equation, and its explicit solution by the Jacobi 
elliptic functions. Using a different idea, we prove an analogous result in higher dimensions for the 
geodesic flow of a closed Riemannian manifold of constant negative curvature. 


1. Introduction 

The goal of this paper is to investigate certain natural but insufficiently explored connections 
between hyperbolic dynamical systems and subriemannian (or Carnot-Caratheodory) geometry. A 
subriemannian geometry on a smooth connected manifold M is a geometry defined by a nowhere 
integrable distribution E, called a horizontal distribution, equipped with a Riemannian metric g. 
Both E and g are required to be at least continuous but in most scenarios they are usually C°°. 
Since we can extend any partially defined Riemannian metric to the entire tangent bundle and the 
extension does not affect the properties of the subriemannian geometry, we will always assume that 
g is defined on the entire tangent bundle. 

If 7 : [a, 6] —)• M is a horizontal (i.e., tangent to E) path, its length is defined in the usual way by 

ItI = [ Il7(i)ll dt, 

J a 

where ||u|| = yg{v^^, for any vector v £ E. 

A horizontal distribution E on M is called nowhere integrable if for every p £ M and every e > 0 
there exists a neighborhood U of p in M such that every point in U can be connected to p by a 
horizontal path of length < e. In particular, every two points of M can be connected by a horizontal 
path. This definition avoids certain undesirable pathological behavior which can arise if E is not 
smooth; see [SimlO]. 

The subriemannian distance between x,y £ M is given by 

dnix^y) = inf{|7| : 7 is a horizontal path from x to y}. 

A subriemannian geodesic from x to p is any horizontal path 7 which minimizes length among all 
horizontal paths connecting x and y. Thus I7I = duix^y). 
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Recall that a C°° horizontal distribution E is called bracket generating if any local smooth frame 
{Xi, ..., Xf;} for E together with all its iterated Lie brackets span the entire tangent bundle of M. 
(In the PDE literature, the bracket generating condition is called the Hormander condition.) By 
the Chow-Rashevskii theorem [Mon 02 ] any bracket generating distribution is nowhere integrable. 

It is sometimes the case that a horizontal distribution E splits into two integrable orthogonal 
distributions, E = Ei (B E2, and E is in turn orthogonal to a globally defined vertical distribution 
V, with TM = E (B V. If E is bracket-generating, then any motion in the vertical direction 
is due to the fact that iterated Lie brackets of vector fields in Ei and those in E2 generate the 
entire tangent bundle. Given a “vertical” curve c tangent to V with endpoints x and y and a unit 
speed subriemannian geodesic 7 connecting x and y, it is natural to ask the following question (see 
Figure 1 ). 


V 



Figure 1. A subriemannian geodesic 7 connecting x = c(0) and y = c(l), where c 
is a path tangent to the vertical bundle V. 


Question. IfTM = E (BV and E = Ei(B E2, how is the energy of a subriemannian geodesic 7 
connecting endpoints of a curve tangent to V shared between its projections onto Ei and E2 ? 


Stated more precisely, if 7 : [ 0 , M is a unit speed horizontal path and F = Fi 0 £"2, then 
f{t) = wi{t) + W2{t), with Wi{t) € Ei. We define 

rl 


£1(1) = IhiWf dt, 


for z = 1 , 2 , and think of Si{'y) as the energy of the projection of 7 to F*. Clearly, 0 < £i{'y) < 1 
and ^1(7) +^2(7) = 1 - If ^^1(7) = ^^2(7) we call 7 a (Fi, F2)-balanced horizontal path. The above 
question therefore asks if every subriemannian geodesic connecting endpoints of a vertical path is 
(Fi, F2)-balanced. 


Example 1 (The Heisenberg group). The Heisenberg group is a subriemannian geometry on M = 
defined by the horizontal distribution F which is the kernel of the 1 -form a = dz — ^{xdy — ydx). 
The Riemannian metric on F is defined by ds'^ = dx"^ + dy^. The vector fields 


y d y d 


d X d 

2 ^ 


form a global orthonormal frame for F. It is not hard to check that [Ai,A2] = ^ =• -^0) so 
F is bracket-generating. Since [Xq,Xi] = [Xq,X2] = 0 , the Heisenberg group is nilpotent. We 
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will show in § 2.2 that every Heisenberg subriemannian geodesic whose endpoints differ only in the 
z-coordinate is (i?i,ill2)-balanced. This follows easily from the fact that Heisenberg geodesics are 
lifts of circles in the xy-plane. 


Subriemannian geometries whose horizontal distributions have a natural splitting into two inte- 
grable distributions occur frequently in hyperbolic and partially hyperbolic dynamical systems. For 
instance, if / is a partially hyperbolic diffeomorphism of a compact manifold M, then / preserves 
two invariant bundles called the stable and the unstable bundles, both uniquely integrable. 
Transverse to them is the center bundle E^, which is not always integrable. Although E^ and E^ are 
usually not smooth, their sum E^'^ = E^^E"^ frequently has the so called accessibility property. This 
means that any two points in M can be joined by a continuous piecewise smooth path whose smooth 
legs are alternately tangent to E^ and E"^. Thus E^^ naturally defines a subriemannian geometry 
on M and we can take H = E^'^ and V = E^. Accessibility plays an important role in partially 
hyperbolic dynamics where it is an essential ingredient in the study of stably ergodic systems and 
the Pugh-Shub conjecture [PS 04 ]. The main difficulty with the subriemannian geometry defined by 
E^^ is that it lacks smoothness, so it is not amenable to analysis using standard techniques. 

In this paper we consider the case of contact Anosov flows, where the natural horizontal distri¬ 
bution is always at least . This is a scenario which is in a sense diametrically opposite to that of 
the Heisenberg group. 

Recall that a non-singular smooth flow = {ft} on a closed (compact and without boundary) 
Riemannian manifold M is called an Anosov flow if there exists an invariant splitting TM = E^^ 0 
EC 0 ^uu that E^ is spanned by the infinitesimal generator X of the flow, E^^ is uniformly 
exponentially contracted and is uniformly exponentially expanded by the flow in positive time. 
We call E^^ and E^"^ the strong stable and strong unstable bundles; E'^ is the center bundle. 

A contact structure on a manifold M of dimension 2 n 0 1 is a hyperplane field E which is 
as far from being integrable as possible [MS 99 ]. This means that there exists a 1 -form a such 
that Ker(a) = E and a A (da)'^ is a volume form for M; a is called a contact form for E. Contact 
structures are always bracket-generating. 

A vector field X is called the Reeb vector field of a if a{X) = 1 and X is in the kernel of da, i.e., 
ixda = 0. An Anosov flow is called contact if FI = E^'^ is a contact structure (and in particular 
C^) and the infinitesimal generator X of the flow is the Reeb vector field for the contact form a 
for E^^ with a{X) = 1 . Our goal is to understand the subriemannian geometry defined by the 
distribution associated with a contact Anosov flow. We will call subriemannian geodesics of 
this geometry sw-subriemannian geodesics. An srt-subriemannian geodesic will be called su-balanced 
if it is balanced with respect to the splitting 0 Fl““. 

Contact Anosov flows have good dynamical properties; in particular, they exhibit exponential 
decay of correlations (cf., [Liv 04 ]). Until recently however, the only known contact Anosov flows 
were the geodesic flows of Riemannian or Finsler manifolds; in [FH 13 ] Foulon and Hasselblatt 
used surgery near a transverse Legendrian knot to construct many new contact Anosov flows on 
3 -manifolds which are not topologically orbit equivalent to any algebraic flow. 

Assume now that $ is a contact Anosov flow on a 3 -manifold M. Denote its infinitesimal generator 
by X and let Y and Z be unit (with respect to some Riemannian metric g whose volume form equals 
the contact volume form) vector fields in E^^ and Fl““ respectively. Then TxftfY) = g{x,t)Y and 
Txft{Z) = \{x,t)Z, for some 1 -cocycles /r, A : M x M ^ M, with A/r = 1 , where Tft denotes the 
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tangent map (i.e., derivative) of the time-t map ft of the flow. Thus 

[X,Y]=aY and [X,Z] = -aZ, (1.1) 

where a{x) = = A(x,0), for all x G M. Since is contact, it follows that [Y,Z] is 

transverse to 

Definition. An Anosov flow ^ on a 3 -dimensional closed manifold will be called a special contact 
Anosov flow if there exists a Riemannian metric g and a global orthonormal frame {X, Y, Z) 
relative to g such that: 

(a) E'^ = MX, E^^ = MT and E^^ = RZ. 

(b) [X, y] = y, [y, Z] = X and [Z, X] = Z. 

Note that X,Y,Z, and g are required to be only C^. The following lemma shows that must in 
fact be C°°. 

Lemma 1 . 1 . If ^ is a special contact Anosov flow, then (with the notation as above), X,Y and Z 
are all C°°. 

Proof. Let be the coframe dual to (X, y, Z). Since X, y and Z are C^, so are a, f and 7 . 

We have: 

1 = a(X) 

= a([y,Z]) 

= ya(Z) - Za(y) - da{Y, Z) 

= -da{Y,Z). 

We can show in a similar way that da{X, Y) = da{X, Z) = 0. Therefore, da is relative to a 
frame, hence a is, in fact, C^. It follows analogously that f and 7 are also C^. Hence X, y and Z 
are all as well. By bootstrap, it follows that X, y and Z are in fact C°°. □ 

Remark. Observe that if $ is a special contact Anosov flow, then X, Y and Z span a copy of the 
Lie algebra s[(2,M), so the universal cover of M is the Lie group SL(2,M) and the lift of the Anosov 
flow to the universal cover is an algebraic one. In particular, M is a quotient of SL(2,M) by a 
discrete cocompact subgroup. In other words, special contact Anosov flows are precisely algebraic 
Anosov flows. 

Ghys [Ghy87] showed that in three dimensions every contact Anosov flow $ with C°° strong 
bundles is (7°° equivalent to an algebraic flow on a quotient X = T \ SL(2,M), in the sense that 
there exists a (7°° diffeomorphism /i : X —>• M that sends the orbits of the “diagonal” flow on 
X to the orbits of <h. Therefore, every contact Anosov flow with C°° strong bundles is C°° orbit 
equivalent to a special Anosov flow. 

Our main results are the following. 

Theorem A. Let $ = {ft} be a special contact Anosov flow on a closed Riemannian 3 -manifold 
M. Then there exists a 6 > 0 such that for all x & M and |t| < 6, every su-subriemannian geodesic 
connecting x and ft{x) is su-balanced. 

In higher dimensions we prove a result analogous to Theorem A if the contact Anosov flow is the 
geodesic flow on the unit tangent bundle of a manifold with constant negative sectional curvature. 
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Theorem B. Let $ = {/^j he the geodesic flow of a closed Riemannian manifold N with constant 
negative sectional curvature on its unit tangent bundle M. Then there exists a 6 > 0 such that for 
all X ^ M and |t| < 6, every su-subriemannian geodesic connecting x and ft{x) is su-balanced. 

Outline of the paper. In Section §2 we review some basic results on Anosov and geodesic flows, 
subriemannian geodesics, and the solution of the harmonic oscillator equation via Jacobi elliptic 
functions. Theorem A is proved in Section §3 and Theorem B in Section §4. We conclude the paper 
with a list of open questions in Section §5. 

Acknowledgments. We are grateful to Alan Weinstein, who generously offered the main idea of 
proof of Theorem B. We would also like to thank him for many inspiring conversations and moral 
support over the years. 


2. Preliminaries 

2.1. Anosov flows. A non-singular smooth flow = {ft} on a closed Riemannian manifold M 
is called an Anosov flow if there exists an invariant splitting TM = 0 such that 

is spanned by the infinitesimal generator of the flow and there exist uniform constants c > 0 , 
0 < /i_ < /i+ < 1 and A+ > A_ > 1 such that for all v € E^^, w € and f > 0, we have 

lli^ll < \\Tft{v)\\ < c/r+ ||u|| , (2.1) 

and 

^A*_ ||'u;|| < \\Tft{w)\\ < cA+ \\w\\ . ( 2 . 2 ) 

The strong stable E^^ and strong unstable bundles are in general only Holder continuous 
[HPS77], but they are nevertheless always uniquely integrable giving rise to the strong stable and 
strong unstable foliations denoted by and respectively. The codimension one distribution 
Esu _ ^ss 0 ^uu jg g 0 nerally not integrable; if it is, then by Plante [Pla72], the flow admits a global 
cross section and is therefore topologically conjugate to a suspension of an Anosov diffeomorphism. 

The bundles = E'^ (B E^^ and are called the center stable and center unstable 

bundles. They are generically only Holder continuous [HPS77], but are always uniquely integrable 
[Ano67]. However, if diniM = 3 and the flow is and preserves the Riemannian volume, then it 
follows from the work of Hurder and Katok [HK90] that and E^“ are both of class and the 
transverse derivatives of both bundles are C^-Hdlder, for all 0 < 9 < 1. 

Without loss we will always assume that all the invariant bundles of an Anosov flow are orientable. 
(If not, we can pass to a double cover.) 

We will need the following easy lemma. 

Lemma 2.1. Let ^ be a contact Anosov flow on a 3-manifold M. Then E““ and E^^ are both . 

Proof. Since E®“ is by assumption and E^^ and E™ are by [HK90], it follows that E^^ = 
E“ n E^“ and E““ = E™ n E^“ are also C^. □ 

Geodesic flows. In this section we briefly review some basic facts about geodesic flows. If N 
is a Riemannian manifold, then its geodesic flow <I> = {fl} restricted to the unit tangent bundle 
M = T^N of A admits a canonical contact form (cf., [Pat99]). If the sectional curvature K oi N 
is negative, then is known to be of Anosov type [Ano67, Ebe73], in which case E®“ is a contact 
structure and TM = E^ © E®“ is an orthogonal splitting with respect to the Sasaki metric [Pat99]. 
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If the sectional curvature K is constant (and negative), then and are C°°, but if K is 
variable, then and are only of class for some 0 < 6* < 1 [HP75]. 

Assume now that K is constant and negative. Without loss we can assume that AT = — 1. Then 
(cf., [Ano67]) there exists a constant c > 0 such that 

\\TU{v)\\=e-^^\\v\\ and \\T h{w)\\ = \\w\\ , 

for all t G M, u G E^^ and w G In other words, the flow contracts all stable directions and 

expands all unstable directions at the same rates at all points of M. For simplicity, we will assume 
that c = 1; this can always be achieved by a constant time change. 

Let E be an isometry of N and denote by iT the restriction of FT to the unit tangent bundle 
M = T^N of N. Then for any unit-speed geodesic t c{t) in N, we have iT(c(t)) = F^,{ftc{0)) = 
/t(F*(c(0)), which implies that iT preserves the geodesic vector field X. We claim that F* also 
preserves the strong stable and strong unstable IT““ foliations of the geodesic flow. Indeed, 
since E^ o ft = ft o F^, for any vi,V 2 in the same IF*^-leaf, we have: 

d{ft{FM),ft{F.iv2))) = d{F,{ftivi)),FMv2))) 

< Kd{ft{vi), ftiv 2 )) 

^0, 

as t —>■ oo, where K is the Lipschitz constant of F* (which is finite, since M is compact and F* is 
smooth) and d denotes the distance function on M induced by the Sasaki metric. Therefore, F*(ui) 
and F*(u 2 ) lie in the same W^^-leaf. Thus F* preserves and TF* preserves F^^. It can similarly 
be shown that F““ is also invariant with respect to TF*. Thus FF* preserves the splitting F'^0F®“. 
An analogous statement is true for any lift F* of F* to the universal Riemannian covering space M 
of M. 

Recall that if FT = —1, then N = H^/T, where T is a group of isometries of H"' acting freely and 
properly discontinuously on it [Boo03] . It is clear that T also acts freely and properly discontinuously 
on the unit tangent bundle F^H” of H"' and that (F^]HI"')/r is isometric to F^(]HI"'/r) = T^N = M. 
Thus M is isometric to the universal Riemannian covering space of (F^BI"')/r. Since ra > 3, F^H” is 
simply connected, so it is the universal covering space of (F^]HI”)/r. Thus M is isometric to F^H”. 

Lemma 2.2. For all u,v € M there exists an isometry F of M such that F(u) = v and F leaves 
the lift X of the geodesic vector field X invariant. 

Proof. By the above observation, M is isometric to F^EI*^, so we can identify u, v with unit tangent 
vectors to H” at some points x,y £ respectively. There exists an isometry / of such that 
f{x) = y and Tf{u) = v (see [Boo03]). Since isometries map geodesics to geodesics, Tf leaves X 
invariant. Thus F = Tf I'T’ien has the desired properties. □ 

2.2. Subriemannian geodesics. In this section we briefly review subriemannian geodesic equa¬ 
tions. We follow [Mon02]. 

Let F be a bracket-generating distribution on a smooth manifold M. For each smooth vector 
field Ai on M we define the momentum function Px '■ T*M —)• M by 

Px{p)=p{X), 

for any p G T*M, where T*M is the cotangent bundle of M. Thus the momentum function of X is 
just the evaluation of any covector on M at A. 
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The subriemannian Hamiltonian H oi E is the map H : T*M —)■ M defined by 

H{p) = ^{p,p), 

where (•, •) denotes the cometric on T*M induced by the Riemannian metric g on E (see [Mon02]). If 
(Xi,..., Xk) is a local horizontal frame and gij = g(Xi,Xj), then the Hamiltonian can be expressed 
as 

H=^2^s‘’Px.Px„ 

hj 

where g^^ are the entries of the inverse of the matrix [gij]- In particular, if {Xi,... ,Xk) is a local 
orthonormal frame for E, then 



2 = 1 

The normal geodesic equation for E is the equation 

f = (2.3) 

where / ; T*M —)• M is a smooth function and {f,E[} denotes the Poisson bracket of / and H. 
Projections of the solutions to (2.3) to M are called normal geodesics. 

Recall that 

{f,H} = cv{Xf,XH), 

where to is the (canonical) symplectic form on T*M and Xf^Xn are the Hamiltonian vector fields 
defined by /, H, respectively. It is well-known that the Poisson bracket defines a Lie algebra structure 
on the ring of smooth functions on T*M and that the map / i— {f,E[} satisfies the Leibniz rule thus 
defining a vector field on T*M (which of course is exactly X//). Recall also that {Pxj ^V} = —P[x,Y]^ 
for any smooth vector helds X, Y on M. 

The equation (2.3) is to be interpreted in the following way: if t i-A p(t) is an integral curve of 
the Hamiltonian vector field Xh and if / : T*M —)• M is any smooth function, then 

= {/.»)&>(*))■ 

In canonical coordinates (xi,..., Xn]Pi, ■ ■ ■ ,Pn), where (xi,..., Xn) are local coordinates on M and 
Pi = PQjQxi, tbe geodesics equations assume the familiar form: 

dH dH 

Xi = Pi = 

opi 

Theorem 2.3 ([Mon02]). Let t T{t) be a solution to the normal geodesic equation (2.3) and let 
7 be its projection to M. Then every sufficiently short arc of j is a subriemannian geodesic. If E 
is a 2-step distribution, then every subriemannian geodesic is normal. 

Recall that H is a 2-step distribution if for any local frame Xi,... ,Xk for E, the vector fields 
Xi,..., Xk together with their first-order Lie brackets [Xi,Xj] (1 < i,j < k) generate the entire 
tangent bundle. 

Assume now dimM = 3 and H is a contact structure. It is easy to see that H is a 2-step 
distribution. Let {Xi,X 2 ) be a local orthonormal frame for E and a a contact form for E. Denote 
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the Reeb field of a by Xq. Clearly, {Xq,Xi,X2) is a local frame for TM. The structure contants 
of the frame {Xq,Xi,X 2 ) are smooth functions defined by 




fc=o 


It follows that 


{PXi,PXj} = - CijPXk ■ 


k=0 

The subriemannian Hamiltonian corresponding to the frame (Xi,X 2 } is 

1 

2 ' 


P = i^iPk+Pk)- 


Introduce fiberwise coordinates {Pxq, Pxi, PX 2 ) on T*M. In these coordinates the normal geodesic 
equations are 


X = PxiXi + PX2X2 
Px, = {Px.,H}, 


for i = 0,1,2. 


Example 2 (The Heisenberg group, continued). We will show that every Heisenberg subriemannian 
geodesic whose endpoints differ only in the ^-component is balanced with respect to the splitting 
E = El (B E 2 . See Example 1. Since {Xi,X 2 ) is an orthonormal frame for E, the subriemannian 
Hamiltonian is 

H = l{Pk+Pk) 

and the subriemannian geodesic equation is / = {/,//}. 

Using [Xi,X 2 ] = d/dz =: Xq, [Xo,Xi] = [Xo,X 2 ] = 0, we obtain {Pxi,Px 2 } = -Pxo, 
{Pxq, Pxi} = {Pxo,Px 2 } = 0. Therefore the subriemannian geodesic equations are 

P = PxiXi + Px 2^2 
Pxo = 0 

Pxi = -PxoPx2 
Px2 = PxqPxi, 

where p = {x,y,z). Since geodesics travel at constant speed, we can restrict the equations to the 
level set P^^ + Px 2 = ^ oi H and reparametrize Pxi and Px 2 by 

Pxi = COS0, Px 2 = sin0. 

It is not hard to check that the last three geodesic equations are equivalent to 

0 = Pxq, Pxq = 0 . 

Thus 0 = 0, so 9{t) = VQt + 6 * 0 , where vq = 0(0) = Pxo(O) and 0o = 0(0)- It follows that every 
Heisenberg geodesic satisfies 

p = cos(vot + do)Xi + sin(uot + vo)X 2 , 

with real parameters vq and 0o as above. Note that x = cos(vot + 0o) and y = sin(uot + 0o). 
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Now assume that a subriemannian geodesic 7 : [0,£] —)• connects two points which differ only 

in the ^-coordinate, i.e., they lie on an orbit of the flow of Xq. Projecting to the xy-plane we obtain 

f cos(vot + 9o) dt = x(i) — x(0) = 0, f sin(uot -|- 9q) dt = y{£) — y(0) = 0. 

Jo Jo 

Thus vq(. must be an integer multiple of 27r. Let us show that 7 is balanced with respect to the 
splitting E = RXi 0 RX 2 . We have 

I ri 1 

Siij) = - / cos'^{vot + 9o) dt and £’ 2 ( 7 ) = 7 / sin^{vot + 9o) dt. 

^ Jo ^ Jo 

Hence 

1 

^ 1 ( 7 ) -S 2 {'y) = J J {cos'^ivot + 9o) - sin^(uot + 0o)} dt 
cos 2{vot + 9o) dt 

= 0 , 

since voi = 27rn, for some integer n. 

2.3. Harmonic oscillator and Jacobi elliptic functions. To make the paper as self-contained 
as possible, we review in some detail the method of explicitly solving the harmonic oscillator (i.e., 
unforced undamped pendulum) equation 

0 + cj2sin6» = O (2.4) 

by the Jacobi elliptic functions sn and cn, defined below. We closely follow [MeyOl], adding results 
we need along the way. 

Let 0 < A; < 1. The Jacobi elliptic functions sn(t, k), cn{t, k) and dn(t, k) are defined as the unique 
solutions x{t),y{t) and z{t) of the system of differential equations 

X = yz 
y = —zx 
z = —k'^xy, 

satisfying the initial conditions 

x( 0 )= 0 , y{0) = l, z( 0 ) = l. 

The paramater k is called the modulus. Some basic properties of sn, cn and dn are listed in the 
following proposition whose proof can be found in [MeyOl] . 

Proposition 1. (a) The Jacobi elliptic functions sn, cn and dn are analytic and defined for all 

real t. 

(b) sn^(t, k) + cifi{t, k) = 1 and k‘^svfi{t, k) + dn^ft, k) = 1 , for allt^M. and 0 < /c < 1 . 

(c) Let K = K{k) > 0 be the unique number such that cn{K,k) = 0 and cn{t,k) > 0, for all 
0 < t < K. That is, K is the time it takes cn{t,k) to decrease to 0 from its initial value 1. 
Then sii{t,K) and dii{t,k) are even about K and cn(t,k) is odd about K. 

(d) sn{t,k) and cn{t,k) are AK-periodic in t and dn{t,k) is 2K-periodic in t. 

(e) The function x{t) = sn{t,k) is the unique solution to the initial value problem 

{xf = (1 — 3:^)(1 — k^x^), x(0) = 0, x(0) = 1. 
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Observe that sn(t,k) and cn{t,k) have the same symmetries with respect to K as sint and cost 
have with respect to 7r/2. 


Corollary 2.4. If 

pb pb 

/ sn(t; k) dt = / cn(t; k) dt = 0, 
J a J a 

then b — a is an integer multiple of 4K. Moreover, 

[ sn(t; k) cn(t; k) dt = 0. 


Proof. The proof follows from parts (c) and (d) of the previous Proposition. The calculations are 
analogous to those proving similar properties for the functions sin and cos. □ 

Now consider the pendulum equation 0 + sin0 = 0, with a; > 0. It is not hard to check that 
the “energy” of the oscillator given by 

1 = +y(1-cos0) 

is constant along solutions. We can rewrite I as 

4 2 

\ , Q 

Let y{t) = sin Then y = —9 cos — = ±0y^l — Squaring both sides and solving for 9“^ from 
the equation for I, we obtain 

(yf = (2.5) 

There are four possibilities. 

Case 1: / = 0. Then 9{t) = 0 (mod 27r) and the pendulum is in the stable downward 
equilibrium. 


Case 2: 0 < / < We look for a solution in the form y{t) = A sn{B{t — to); k), for some 
constants A, B, to and 0 < A; < 1, and obtain 

y{t) = k sn{uj{t — to);k), with fe = 

OJ 

Therefore, for any to G K, 

9{t) = 2 arcsin {k sn(cu(t — to); k)} , 

is a solution to (2.4). This case corresponds to the pendulum swinging back and forth. 
Given a particular solution 9{t), we can compute to using 9(to) = 0. 

Case 3: I = uj^. If 9{0) = 0, then 9{t) = tt (mod 27r) and the pendulum is in the unstable 
upward equilibrium. If 0(0) / 0, then 

y = ±a;(I -y^), 

and it is not hard to check that 


y{t) = ±tanh(a;(t — to)) 

satisfies (2.5) for any to, so the solution to the pendulum equation in this case is 

9{t) = ±2 arcsin tanh(ti;(t — to)). 
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This solution is a saddle connection connecting two downward equilibria. 

Case 4: I > We seek a solution in the form y(t) = sn(74(t — to)] k), for some constants 
A, to and 0 < A; < 1, and obtain 


y{t) = sn{'/l{t — to)]k), where k =—j=, 

\/I 


which yields 


9{t) = 2 arcsin ^sn(\/7(t — to); k)'j , 

for any to G M. Given a particular solution 9{t), we can compute to from the equation 
0{to) = 0. This case corresponds to circulating orbits, where the pendulum has enough 
energy to go over the top. Note that we need to keep changing the branches of arcsin in 
order to keep 9 increasing or decreasing. 

Lemma 2.5. Let 1 1 —)• 9{t) be a solution to 9 + sin0 = 0, for 0 <t < i. Suppose that 


9{t) , 9{t) , 

/ sin dt= cos — 7 ^ dt = 0. 

Jo 2 Jo 


10 ^ 70 2 

Then: 

(a) I > uj^. 

(b) i^/l is an integer multiple of the period of sn( • ; k), with k = u/y/l. 

(c) Furthermore, 


( 2 . 6 ) 


/ 

70 


sin0(t) dt = 0. 


Proof. The assumption (2.6) implies that only Cases 2 and 4 above are possible. If 

sin —^ = k sn{u}{t — to)] k) 

as in Case 2, then by Proposition 1(b), |sin0(t)/2| < k, so |cos0(t)/2| > y/l — kf, contradicting 
the second part of (2.6). Therefore, only Case 4 is possible, which implies I > a;^, proving (a). 
Furthermore, we know that 

9{t) 


Sin • 


= sn 


I{t - to)]k'j , 


where k = ujj\fl. Thus 


cos ■ 


9{t) 


= ± 1/1 — sin" 


m 


= — sn^ i^yi{t — to); 

= ±cn ^\/7(t — to)] A:'^ 

/ sm- - dt = —= / sn(t; k) dt 

Jo 2 VlJ~V7to 

ri I |>vT(£-to) 

/ COS dt = ±—p / cn(t; k) dt, 

Jo ^ yl J—y/l tQ 

Corollary 2.4 and (2.6) imply that iy/l is an integer multiple of 4iF, where AK is the period of 
sn( • ; k) and cn( ■ ]k). This proves (b). 


Since 


and 
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Finally, again using Corollary 1, we obtain (c): 

pi pi 

/ suiO{t) dt = ±.2 / sn{Vl{t — to);k) cn{Vl{t — tQ)-,k) dt 
Jo Jo 

2 i-Viii-to) 

= ±—p / sn(t; fc) cn(t; fc) dt 

V k J — y/l to 

= 0 . □ 

3. Proof of Theorem A 

Let be a special contact Anosov flow on a closed 3-manifold M. Then 

[X,Y]=Y, [Y,Z]=X, and \Z,X\ = Z, (3.1) 

where Y £ Z £ are (7°° vector fields and (A, T, Z) is an orthonormal frame with respect to 
a fixed Riemannian metric on M. On each fiber of T*M we introduce the coordinates {Px, Py,Pz), 
where Px, Py, Pz are the momentum functions of X, Y, Z, respectively. The subriemannian Hamil¬ 
tonian is 

H = i(p2+p|). 

The subriemannian geodesic equations in these coordinates are 

X = PyYx + PzZx 

Px = {Px,P} 

Py = {Py,P} 

Pz = {Pz,H}. 

Using (3.1) it is not hard to see that 

{Px,H} = Pl-P^, {Py,H} = PxPz, and {Pz, H} =-PxPy- 
Thus the geodesic equations are 

X = PyYx + PzZx 

Px = Pi- Pi 
Py = PxPz 
Pz = -PxPy- 

Since H is an integral of motion, we will consider only solutions lying on the level set ^ defined by 
Py -|- P| = 1. Each fiber of .if is topologically a cylinder x M. We set 

Py = cos 9, Pz = sind, 

so that on .if we have local coordinates (x;d,Pv), where x £ M. In these coordinates the subrie¬ 
mannian Hamiltonian equations become 

X = (cos 9)Y + (sin 9)Z 

0 = -Px 

Px = sin^ 9 — cos^ 9. 


Therefore, 9 satisfies 


9 — cos 29 = 0. 
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Let X £ M be arbitrary and pick a flowbox U for X containing x. Let r > 0 be small enough so 
that y = frix) lies in U and let 7 be a unit speed subriemannian geodesic such that 7 ( 0 ) = x and 
7 (£) = y, where I = dsu{x,y). Then 


7(t) = cos6l(t)y^(t) + sin6»(t)Z^(t), 
lillator equation 6 — cos 29 = 0. 1 

pi pi 

/ cos^ 9{t) dt= sin^ 9{t) dt, 
Jo Jo 


where 9 satisfies the harmonic oscillator equation 9 — cos 29 = 0. Observe that 7 is su-balanced if 

r£ r£ 


which is equivalent to 


f cos 20 (t) 

Jo 


dt = 0 . 


TT 


The substitution cj) = 29 — — in the 0-equation converts it to the standard form 

(f) + 2sin(f) = 0, 

with CJ = y/2. 

Now consider the orbit space S = [//~ of the Anosov flow in [/, where z ~ u; if u; is on the orbit 
of z. It is clear that S is diffeomorphic to an open set in Denote by 7 by the natural projection 
of 7 to S and by Y, and Z the projections of Y and Z along 7 to S. Then relative to the frame 
{Y, Z) along 7 , we have 

7 (t) = (cos 0 (t), sin 0 (t)), 
for 0 < t < £. Since ^{l) = 7 ( 0 ), it follows that 

f cos9{t)dt= f sin9{t) dt = 0. 

Jo Jo 


1 TT 

Using 9{t) = substituting, we obtain 


(j){t) 


f Mt) r 
/ cos dt= sin dt = 0 . 

Jo 2 Jq 


(3.2) 


Lemma 2.5 yields 


fi 

/ sin(/)(t) 

Jo 


dt = 0 . 


But sin (j) = — cos 29 = sin^ 9 — cos^ 0, so 


pi pi 

/ sin^ 9{t) dt= cos^ 9{t) dt = 0, 

JO Jo 

proving that 7 is an su-balanced geodesic. 


4. Proof of Theorem B 

Assume <I> = {ft} is the geodesic flow of a closed Riemannian manifold N with constant negative 
curvature; <I> is defined on M, the unit tangent bundle of N. Without loss we can assume that 

ll^/i(^^)ll = e"* ||u|| and ||r/i(u;)|| = e* ||u;|| , 
for all V £ w £ and all real t. 

Let p & M and r > 0, and set q = fr{p)- Let 7 : [0,.^] —>• M be a unit speed stt-subriemannian 
geodesic from p to g so that I 7 I = i = dguip, q)- We assume r is sufficiently small so that the closed 
curve consisting of 7 and the arc of A-orbit {ft{p) : 0 < t < r} is homotopically trivial. Clearly, 
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there exists a uniform <5 > 0 such that this is true for every two points on the same $-orbit which 
are ( 5 -close (with respect to the Riemannian distance) and for every su-subriemannian geodesic 
connecting them. 

We claim that 

dsuiftip), ft{q)) = dsu{p, q), ( 4 . 1 ) 

for all t G M. 

To prove this, first denote by p, q the lifts of p, q lying in the same fundamental domain of the 
universal Riemannian covering tt : M — >• M. Let t G M be arbitrary and let x = ft{p) and y = ft{q)- 
Denote by x,y the corresponding lifts to the universal covering, so that x = ft{p) and y = ft{q)- 

By Lemma 2.2 there exists an isometry F oi M such that F{p) = x and TF{X) = X, where X 
denotes the lift of X to M. Denote by 7 the lift of 7. Then T o 7 connects x and y and has the 
same length as 7 as well as 7. 




Figure 2 . Proof of Theorem B. 

We claim that its projection 7 = 710^07 is an su-geodesic connecting x and y. See Figure 2 . 
First observe that since F preserves X, and and X are orthogonal, 7 is an su-path. Suppose 
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that it is not minimal. Then there exists a horizontal path c joining x and y such that |c| < |7|. 
This path lifts to a path cm. M connecting x and y, and tangent to the lift of c has the same 
length as c. It follows that vr o F~^ o c is a horizontal path joining p and q whose length is 

o F~^ o c| = |c| < I7I = I7I , 


contradicting the assumption that 7 is a geodesic. This proves ( 4 . 1 ). 

Now let us show that 7 is su-balanced. Write 7(t) = Ws{t) + Wu{t), with Ws{t) £ and 
Wu{t) € E^'^, for all 0 < t < £. Then: 


A. 

dr 


l/r0 7| 

r=0 


A. 

dr 

A_ 

dr 

A_ 

dr 

A_ 

dr 


r=0 


r=0 


^=0 *^0 


^=0 ^0 


[ \\Tfr{ 7 {t))\\ dt 

Jo 

[ \\Tfriwsit)) + Tfr{Wuit))\\ dt 

Jo 

\J\\Tfr{Wsit))f + \\Tfr{Wu{t))f dt 

[ g2r ||u;^(i)||2 (If 

Jo 


2 j |l|u^«(0ll^ - ll^^s(^)ll^} dt 

2^{T„(7)-T,(7)}. 


( 4 . 2 ) 


Here we used the fact that 7 is unit-speed and Ws{t) and Wu{t) are orthogonal. If 7 is not su- 
balanced, i.e., ^^(7) / then by ( 4 . 2 ), there exists t / 0 such that |/t0 7| < I7I, which 

yields 

dsuiftip), Mq)) < Iftoil < I 7 I =dsu{p,q), 


contradicting ( 4 . 1 ). This completes the proof. 


5. Open questions 

Virtually any question involving subriemannian geometry defined by the stable and unstable 
bundles of a partially hyperbolic dynamical system or an Anosov flow is open. Here we list only a 
few. 

( 1 ) If $ is an arbitrary contact Anosov flow on a 3 -manifold, what can be said about £3(7) and 
Sui^), for an su-subriemannian geodesic 7 whose endpoints lie on an orbit of the flow? Are 
there lower and upper bounds for Eg (7) /(7) independent of 7? The same question can be 
asked about arbitrary contact Anosov flows in any odd dimension. 

( 2 ) What if is the geodesic flow of a closed Riemannian manifold with variable negative 
sectional curvature? 

( 3 ) Suppose is a transitve Anosov flow such that has the accessibility property but is 
not necessarily contact. Do su-subriemannian geodesics exist? If so, how close are they to 
being sn-balanced (when their endpoints lie on the same orbit of <h)? 

( 4 ) The same question can be asked for an accessible partially hyperbolic diffeomorphism. 
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